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We consider the Coulomb correction (CC) to the e + e~ pair production related to multiphoton 
exchange of the produced with nuclei. The contribution of CC to the energy distribution of e + 
and e~ as well as to the total pair production cross section are calculated with an accuracy of the 
order of 1 %. The found correction to the total Born cross section is negative and equal to —25 % 
at the RHIC for Au-Au and -14% at the LHC for Pb-Pb collisions. 



Introduction. Two new large colliders with relativis- 
tic heavy nuclei, the RHIC and the LHC, are scheduled 
to be in operation in the nearest future. The charge num- 
bers Z\ = Z2 = Z of the nuclei with masses Mi = M2 = 
M and their Lorentz factors 71 = 72 = 7 = E/M are the 
following 

Z = 79 , 7 = 108 for RHIC (Au - -Au collisions) 

Z = 82 , 7 = 3000 for LHC (Pb - -Pb collisions) . (1) 

Here E is the heavy ion energy in the c.m.s. One of the 
important processes at these colliders is 



Z X Z 2 — » ZiZ 2 e + e 



(2) 



Its cross section is huge. In the Born approximation (see 
Fig. [I] with n = n' = 1) the total cross section according 
to the Racah formula jjj is equal to o~ Bom = 36 kbarn 
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FIG. 1. The amplitude M nn i of the process (2) with n (n) 
virtual photon emitted by the hrst (second) nucleus. 



for the RHIC and 227 kbarn for the LHC. Therefore it 
will contribute as a serious background to a number of 
experiments, besides, this process is the leading beam 
loss mechanism (for details see review Q). 



The cross sections of the process (g) in the Born ap- 
proximation are known with accuracy ~ I/7 2 (see, for ex- 
ample, Rcfs. jl],H and more recent calculations reviewed 
in Refs. HQ). However, besides of the Born amplitude 
AfBom = M11, also other amplitudes M nn i (see Fig. |l|) 
have to be taken into account for heavy nuclei since in 
this case the parameter of the perturbation series Za is 
of the order of unity. Therefore, the whole series in Za 
has to be summed to obtain the cross section with suffi- 
cient accuracy. Following Ref. J5[, we call the Coulomb 
correction (CC) the difference Arcoui between the whole 
sum da and the Born approximation 



da, 



Coul • 



(3) 



Such kind of CC is well known in the photoproduc- 
tion of e + e~ pairs on atoms (see Ref. M and §98 of 0). 
The Coulomb correction to the total cross section of that 
process decreases the Born contribution by about 10 % 
for a Pb target. For the pair production of reaction (0) 
with Z\a <C 1 and Z2Q! ~ 1 CC has been obtained in 
Refs. j?],U. Recently this correction has been calculated 
for the pair production in the collisions of muons with 
heavy nuclei The results of Refs. agree with 

each other in the corresponding kinematic regions and 
noticeably change the Born cross sections. Formulae for 
CC for two heavy ions were suggested ad hoc in Sect. 7.3 
of . However, our calculations presented here do show 
that this suggestion is incorrect. 

In the present paper we calculate the Coulomb correc- 
tion for process (J2j> omitting terms of the order of 1 % 
compared with the main term given by the Born cross 
section. We find that these corrections are negative and 
quite important: 

CTCoui/ffBorn = -25% for RHIC, 



1 



o 



-14% for LHC. 



(4) 



This means that at the RHIC the background process 
with the largest cross section will have a production rate 
25 % smaller than expected. 

Our main notations are given in Eq. (|l]) and Fig. [j], 

besides, (Pa + P 2 ) 2 = 4P 2 = 4 7 2 M 2 , q % = q*) = 
Pi — P-, e — e + + E- and 



CO 



aAZ 2 Z 2 T , PiP 2 



L = In- 



2MiM 2 



ln 7 2 



(5) 



where m is the electron mass. The quantities qjj_ and 
p±± denote the transverse part of the corresponding 
three-momenta. Throughout the paper we use the well 
known function || 



f( Z ) = z w y: 



^ n(n 2 + Z 2 a 2 ) ' 

n— 1 v / 



(6) 



its particular values for the colliders under discussion are 
/(79) = 0.313 and /(82) = 0.332. 

Selection of the leading diagrams and the structure of 
the amplitude. Let M. be the sum of the amplitudes 
M nn i of Fig. [l| It can be presented in the form 



M = Yj M nn> = M B orn + Ml + Ml + M 2 , (7) 

nn'>l 

Ah = Mm' , Mi = M m . M 2 = ^2 M nn> ■ 

n'>2 n>2 nn'>2 

The Born amplitude Mbo™ contains the one-photon 
exchange both with the first and the second nucleus, 
whereas the amplitude Mi (Mi ) contains the one-photon 
exchange only with the upper (lower) nucleus. In the last 
amplitude M 2 we have no one-photon exchange. Accord- 
ing to this classification we write the total cross section 
as 



V\ + 01 + CT 2 



(8) 



where 



oc |Mb orn | 5 
d^ cx 2Re(M Bo rnM*; 
d&! cx 2Re(M B or„M*; 



|M!| 2 . 
|M!| 2 . 



da 2 oc 2Re (M Bo r„M 2 * + M^* + MiM 2 * 



+MiM 2 *J + |M 2 | 2 . 

It is not difficult to show that the ratio o~i/o~Bom is 
a function of {Za) 2 only but not of Za itself. Ad- 
ditionally we estimate the leading logarithms appear- 
ing in the cross sections <7j. The integration over the 
transfered momentum squared q\ and q 2 results in two 
large Weizsacker-Williams (WW) logarithms ~ L 2 for 



the CTBom, m one large WW logarithm ~ L for o\ and 
o i . The cross section cr 2 contains no large WW loga- 
rithm. Therefore, the relative contribution of the cross 
sections <7j is o~\ jo-Barn = o'l/o'Born ~ {Za) IL and 
CT 2 /cTBor„ ~ (Za) 2 /L 2 < 0.4 % for the colliders (@). As a 
result, with an accuracy of the order of 1% we can neglect 
cr 2 in the total cross section and use the equation 



CBo 



(9) 



With that accuracy it is sufficient to calculate o\ and 
g i in the leading logarithmic approximation (LLA) only 
since the next to leading log terms are of the order of 
(Za/L) 2 . This fact greatly simplifies the calculations. 

The calculation in the LLA can be performed us- 
ing the equivalent photon or WW approximation. The 
main contribution to o\ and a± is given by the region 
(wi/7) 2 < -q\ < to 2 and (w 2 / 7 ) 2 < -qt < m 2 , 
respectively. In the first region the main contribution 
arises from the amplitudes Mb hi + Mi (in the second 
region A/Bom + Mi). The virtual photon with four- 
momentum qi is almost real and the amplitude can be 
expressed via the amplitude M 7 for the real photoproduc- 
tion 7 Z 2 — > Z 2 e + e~ (see, for example, §99 of Ref. ||) 



Mi 



Born 



4 ;;aZiiH4— M 



(-if) Wl 



7 • 



(10) 



The amplitude M 7 has been calculated in Ref. ||. We 
use the convenient form of that amplitude derived in the 
works § and 0: 



M 1 = (A M 7 



i/ 2 AM 7 )e i4 



(11) 



where M®° rn is the Born amplitude for the jZ 2 — > 
Z 2 e + e~ process. This Born amplitude depends on the 
transverse momenta p±± only via the two combina- 
tions A = £ + — and B = £+p+± + £-P-J_ where 
£ ± = m 2 /(m 2 + Pj-j_)- The quantity AM 7 is ob- 
tained from M? orn replacing A — ► £ + + — 1 and 

B -> £+P+_L -C-P--L- 

All the nontrivial dependence on the parameter Z 2 a = 
v are accumulated in the Bethe-Maximon phase 



<E> = v In 



and in the two functions (with z = 1 — (— <z 2 /?ri 2 )£+£ 

1 - x 



(12) 



= F(ii/,-h/; l;z) 



(13) 



The function fi(z) and its derivative /{(z) are given 
with the help of the Gauss hypergeometric function 
F(a, b; c; z). 

It can be clearly seen that in the region p±j_ ~ m 2 
the amplitude M 7 differs considerably from the Mr orn 
amplitude and, therefore, the whole amplitude M differs 



2 



from its Born limit Msorn- Let us stress that just this 
transverse momentum region p^_^ ~ m 2 gives the main 
contribution into the total Born cross section CTBom and 
into o~ i . 

Outside this region the CC vanishes. Indeed, for 
p± ± <C m 2 or p^_ ± 3> m 2 the variable z « 1, therefore, 
/i w 1, / 2 w and 



M Bot „ + Mi = M B C 



(14) 



Note that the region p± ± m2 gives a negligible contri- 
bution to the total cross section a, however, this region 
might be of interest for some experiments. 

The results of Ref. || which are used here in the form of 
Eqs. (|l(|)-(|l4|) are the basis for our consideration. These 
results were confirmed in a number of papers (see, for 
example, Refs. |l^ , [l0[|) using various approaches. 

Recently in Refs. JlJ the Coulomb effects were studied 
within the frame-work of a light-cone or an eikonal ap- 
proach. However, the approximations used in Refs. fTl| l 
fail to reproduce the classical results of Bethe and Max- 
imon ||. To show this explicitly, we consider the sim- 
ple case Z\a -C 1, Z 2 a = v ~ 1 in which the princi- 
pal result of Refs. jll) for the amplitude takes the form 
M = M Bo rn + Mi = A/ Bo rnexp(i^) with # = ^lnq 2 .^ 
in obvious contradiction to Eqs. (^0|)-(|ll|). Since in the 
works flTl] | different statements on the applicability range 
of their results can be found, we take as an example the 
common region (uji/ji) 2 <C q 2 j_ <C m 2 . But even in that 
region their expression for the matrix element does not 
reproduce Eq. ( |l4"| ) since their phase \& does not coincide 
with the Bethe-Maximon phase $, i.e. W ^ <f). 

CC to the energy distribution and to the total cross 
section. As it was explained in the previous section, the 
basic expression for the cross section dai in the LLA can 
be directly obtained using the WW approximation. To 
show clearly the terms omitted in the LLA, we start with 
a more exact expression for do\ derived for the case of ixZ 
collisions considered in Ref. || . The reason is that for the 
most interesting region (when the energy of relativistic 
e pairs is much smaller than the nucleus energy) the 
muon in the ijlZ scattering as well as the upper nucleus 
of the ion-ion collision can be equally well treated as 
spinlcss and pointlike particles. 

Using Eqs. (14) and (17) from Ref. || (given in the 
lab frame of the muon projectile on a nucleus target) 
and the invariant variables x± — (j)±P 2 ) / {qiP 2 ) , V — 
{Q1P2)/ '{Pi Pi) w e obtain do~\ for the pair production in 
Z\Z 2 collisions in the invariant form (and at y <C 1) 



do\ 



ia /(Z 2 ){[(l+£) a -l]lni±^ 



4 — a 1 dy 



-dx+dx-8(x + + X- — 1) (15) 



with a = 2(1 + x\ + x 2 _) , £ = (Miy/mf x+X- 
main contribution to o~\ is given by the region 



The 



MfMj 

(W 



1 



(16) 



The corresponding expression for d&\ can be obtained by 
making the replacements 



do' 1 = do-\{q\ — > q 2 , P\ <-> Pi, Z Y <-» Z 2 



(17) 



Below we consider only the experimentally most inter- 
esting case when in the collider system (71 = E\/M\ ~ 
72 = E 2 /M 2 ) both e + and e~ are ultrarelativistic (e± ^> 
m). We assume that the z-axis is directed along the ini- 
tial three-momentum of the first nucleus Pi. 

To obtain the energy distribution of e + and e~ in the 
LLA we have to take into account two regions p± z S> m 
and (— p± z ) S> m where the lepton pair is produced either 
in forward or backward direction. In the first region we 
have x± — e±/e, y — e/Ei, and from Eq. (|l~5|)-(|l6|) we 
obtain in the LLA 



doV = 



(to 7 i) 2 de + de_ 
In - — , (18) 



e + £- 



m <C e± <C m7i . 



In the second region we have x± 
m 2 £/{AEi£ + S-)) and 



da\ 



(2)_ 



4a f(Z 2 ) 1 



4e + e_ 

"Sk 2 " 



In 



7 2 e + £r_ de + de. 



(19) 



m -C e± -C m72 . 
Summing up these two contributions, we find 

4e + e 



da 1 = -8a f{Z 2 ) 1- 



3, 2 



ln-^^. (20) 



To obtain u\ we have to integrate the expressions 
and (|l^) over e_ (with logarithmic accuracy) 



do\ 



28 , (m 7 i) 2 de+ 
—o-oJ{Z 2 ) In g ~ ' 



9 

m <C e + -C JT171 , 



a (2) 28 tt v s 1 (7i£+) 2 

o^i = — 7rcr f{Z 2 ) In - 

9 m z e + 

m -C e+ -C m-72 

from which it follows that 

28 2 cfc + 
d(Ti = — —o- j{Z 2 ) M71 . 

y E-\. 



(21) 



(22) 



(23) 



The further integration of Eqs. (|2l|), ( |22|) over e + results 
in 



0-1 



28 



o-of{Z 2 



In 



P1P2 
2M 1 M 2 



(24) 
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This expression is in agreement with the similar result 
for the /iZ scattering (see Eq. (31) from S for Z\ = 
l,Z 2 = Z). 

The corresponding formulae for d\ can be obtained 
from Eqs. ©, © and (||) by replacing 71 <-> 72, Zi 
Z2. The whole CC contribution <i<7coui = d(<7i + 5i) for 
the symmetric case Z\ = Zi = Z and 71 = 72 = 7 takes 
the following form 



oul 



-16<7„/(Z) 1 



4e + e_ 



de+ds- 



at m C £± < 7717 , 

docoui = -- 
at m -C e+ <C TO7 , and 



56 



CTCoul = —-g-<7of(Z) L 2 



(25) 



(26) 



(27) 



The size of this correction for the two colliders was 
given before in Eq. ([!]) . The total cross section with and 
without Coulomb correction as function of the Lorentz 
factor 7 is illustrated in Fig. @ for Pb nuclei. 



" ^Born 

" ^Born + Ocoulomb 




50 100 200 500 1 3 2- 1 3 5-1 3 
7 = E/M 

FIG. 2. The total cross section of the process ZZ — » 
ZZe + e~ with (solid line) and without (dashed line) Coulomb 
correction as function of the Lorentz factor 7 of Pb nuclei 
(Z = 82). 



ourselves to the Coulomb corrections for the energy dis- 
tribution of electrons and positrons and for the total cross 
section. In our analysis we neglected contributions which 
are of the relative order of ~ (Za) 2 /L 2 . The CC to the 
angular distribution of e + e~ can be easily obtained in a 
similar way, however only with an accuracy Za/L 2 . 

Since our basic formulae ([l5|) , ( |l7j ) are given in the in- 
variant form, a similar calculation can be easily repeated 
for fixed-target experiments. This interesting question 
will be considered in a future work. 
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Conclusion. We have calculated the Coulomb correc- 
tions to e + e~ pair production in relativistic heavy ion 
collisions for the case of colliding beams. Our main re- 
sults are given in Eqs. (p5|)-(p7|). We have restricted 
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